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Abstract 

We construct the low-energy effective theory for the SO (5) model of high-T c super- 
conductivity, recently proposed by S.C. Zhang ( |cond-mat /96 1 140|) . This permits us to 
develop a systematic expansion for low-energy observables in powers of the small symmetry- 
breaking interactions. The approximate SO (5) symmetry predicts relations amongst these 
observables, which are model-independent consequences of Zhang's proposed symmetry- 
breaking pattern. 
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1. Introduction 



In a recent remarkable paper [1] Zhang has suggested an elegant framework for under- 
standing the relationship between the superconducting (SC) and antiferromagnetic (AF) 
properties of the high-T c superconductors. His proposal comes in two parts. First, he 
argues that these systems enjoy an approximate SO (5) symmetry which contains as a 
subgroup the SO (3) x SO (2) symmetry of spin rotations and electromagnetic gauge trans- 
formations. The SO (5) symmetry is only approximate, in the sense that it is explicitly 
broken to SO (3) x SO (2) by small interaction terms in the Hamiltonian, whose scale we 
denote by E^. It is also explicitly broken by the doping of electrons or holes away from 
half filling. Second, the 50(5) symmetry is argued to be spontaneously broken to 50(4) 
by the dynamics which binds the electrons into spin-singlet pairs. The energy scale, E 5: 
which characterizes the order parameter for this breaking is imagined to be comparable to 
the Fermi energy, and is taken much larger than E s ^. For energies much smaller than E§, 
the symmetry-breaking dynamics is purely concerned with how the 5'0(5)-breaking order 
parameter aligns relative to the direction of explicit symmetry breaking. In ref. [1] Zhang 
argues convincingly how such a framework can synthesize a great many features of these 
systems, and proposes a model low-energy effective theory which qualitatively describes 
their properties, and does so quantitatively in the vicinity of the critical point between the 
antiferromagnetic (AF) and superconducting (SC) regimes. 

Here we generalize Zhang's model to incorporate the most general possible interactions 
which are consistent with the assumed symmetry-breaking pattern. Being the most general 
possible such lagrangian, it must then incorporate the low-energy limit of any particular 
microscopic model which realizes these symmetries. Because the low-energy behaviour 
involves the interactions of Goldstone (and pseudo-Goldstone 1 [2]) bosons, they are guar- 
anteed to interact only weakly at low energies. Consequently, a mean-field treatment of 
the effective theory is justified to describe these systems for energies well below E 5 , and 
away from any critical points. This description closely parallels that of chiral perturbation 
theory [3], which has been very successfully used to describe the low-energy interactions 

1 A pseudo-Goldstone boson is the Goldstone boson for a symmetry which is only approximate. 
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of pions and nucleons within the framework of Quantum Chromodynamics. 

Eq. (5) expresses our main results, which consist of a number of relations amongst the 
low-energy observables of these systems which are generic consequences of, and therefore 
strong tests of, the approximate SO (5) symmetry. For the sake of brevity we here simply 
state our results, with a more thorough discussion to appear elsewhere. 

We start with the lagrangian density for the four would-be Goldstone bosons, 6> Q , 
a = 1,...,4, corresponding to the spontaneous breakdown SO (5) — > 50(4) [5]. Writing 
this as a derivative expansion, and, for simplicity, ignoring the system's spatial anisotropy 
gives: 2 

f 2 f 2 

£inv(#) = -V + — d t n T d t n - J -f Vn T • Vn + higher derivatives. (1) 

Here Vo, ft and f s are real constants which are of order unity (in units where E§ = 1). 
The four variables, 9 a , parameterize the ground-state field configurations, which fill out 
the space 50(5)/ 'SO (4) or, equivalently, the 4-sphere. We use in (1) coordinates consisting 
of a five-dimensional vector — the 'superspin' — n(9), having unit length, n T n = 1. n 
transforms simply with respect to 50(5): n — > gn, for g a real, orthogonal five-by-five 
matrix. 

In order to explicitly break 50(5) down to 50(3) x 50(2) we introduce the symmetry- 
breaking matrix: M = ediag (m q , m s , m s , m s , m q ). Here e = E^/E^ ~ few %, is the 
small dimensionless parameter which describes the strength of the explicit 50(5) breaking. 
Following Zhang, we choose the spin 50(3) to rotate ri2, and into one another, while 
the electromagnetic 50(2) rotates ri\ into n 5 . The explicit breaking of 50(5) symmetry 
by doping is achieved within the microscopic theory by coupling a chemical potential, 
/U, to the electron's electric charge, Q. In units of E 5 , we expect [i to lie in the range 
\[x\ < 0.1 — 0.2 for the doping of high-T c systems. Since electric charge is a generator of 
50(5), we take: Qij = q(Su — S±j 5$i). q = ±2 is the electric charge (in units of e) of 
the superconducting order parameter. 

The lagrangian for the system is then the most general local function of n(6>), M and 

We choose units for which h=c=k B =E 5 =l. 
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fiQ which satisfies C(gn, gMg T , gnQg T ) = C(n, M, nQ). The leading terms in a derivative 
expansion therefore are: 



= -V + f t [Ad t n Q T d t n Q + B d t n s T d t n s + C (n Q T d t n Q ) 2 ] (2) 
- fl DVn Q T ■ Vn Q - EVn s T ■ Vn s - F (n Q T Vn Q ) 2 j + higher derivatives, 

where n Q and n s respectively are the components of n within the electromagnetic and spin 
subspaces. The functions V, A, B, C, D, E and F are all functions of n, M and \xQ. For 
instance the scalar potential, V, has the general form: 

oo 

V[n, M, nQ] = 2" fc (vfc Q n Q T n Q + v ks n s T n s ) k . (3) 
k=i 

Analogous expansions may also be written for the functions A through F (with coefficients 
denoted cik Q ,aks through fkQifks)- Coupling to electromagnetic fields, A^, is achieved 
by replacing ordinary derivatives of n Q with covariant ones: d^n Q — > D^n Q = d^n Q — 
qeA^Qrig. 

It is useful to express the conserved spin and electromagnetic currents, which are 
obtained from (1) and (2) using Noether's theorem, for the effective theory, since it is to 
these that external probes often couple. The terms involving the fewest derivatives are: 

Pem = -ft (l + 2A)n Q T d t n Q , j em = f 2 s (l + 2D)n Q T QVn Q , 

V ' V J (4) 

Pspin = ft (l + 2B)n s x dtn s , Jspin = - f s (l + 2E)n s x Vn s . 

Notice that eqs. (1), (2) and (4) also capture the form for the free energy and currents 
at nonzero temperature if all coefficients are understood as functions of T. This is because 
these expressions are the most general possible consistent with the derivative expansion 
and the assumed symmetry-breaking pattern. 

The predictive power of eq. (2) emerges once we perturb in the small symmetry- 
breaking parameters e and fx. This requires knowing the dependence of the coefficients 



in eq. (3) (and its analogs for the functions A to F) on these small parameters. It is an 
easy exercise to see that the following properties hold: Vk S , Uksibks, Cks, dks, e ks an d fks 
are independent of \i. These coefficients are at most 0(e) in size, as are all of the others 
when restricted to half-filling: fx = 0. The antisymmetry of the matrix Q further implies 
that vuq, afc Q , 6fc Q , dk Q and ekQ must be even under the interchange fx — > — /z. The same 
is not true for the coefficients Ck Q and fk Q . Because most low-energy observables do not 
appreciably depend on these last coefficients, they also enjoy a symmetry between positive 
and negative doping. 

Perturbing in e and fi is powerful because all interactions amongst the 9 a in eqs. (1) 
and (2) vanish in the limit e = fx = 0. Mean field theory is therefore justified so long as 
three criteria hold: (i) fx and e are small; (ii) the derivative expansion is justified (long 
wavelengths compared to 1/E$; and (Hi) we are not in the vicinity of critical points, for 
which non-Goldstone modes become massless, with the associated uncontrollable fluctua- 
tions in the infrared. 

We now summarize some of the predictions which follow to leading order in e and \x. 
For details concerning their derivation, and some results beyond leading order, see [4]. To 
leading order it suffices to take A = B = C = D = E = F = and keep only the 50(5)- 
invariant interactions having two derivatives. The symmetry-breaking terms then enter 
predictions only through the potential, V. It suffices to keep only the terms having k = 1, 2 
in the expansion (3), and to expand the coefficients v\ Q and V2 Q to linear order in fx 2 : e.g.: 
viq = i>i Q + vl Q fx 2 etc.. In this limit the phase diagram, long- wavelength dispersion 
relations, and electric and magnetic stiffness are all determined by five parameters, which 
we choose to be m 2 = v^ Q - v ls = 0(e), k = -v\ Q = 0(1), f = (v^ Q ) 2 = 0(1), f s = 0(1) 
and f t = 0(1). 

Since these five parameters relate more than five low energy observables for the SC 
and AF phases, we obtain testable relationships amongst these observables, which we now 
summarize. 

Writing the magnitudes of n s and n Q as \n s \ = sin 9 and \n Q \ = cos 6, and minimizing 
the potential with respect to 9 gives two generic phases: 6> m i n = 0, tc (SC); and 6> m i n = f, 
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(AF). The curvature of the potential at these extrema are M 2 F = (d 2 V / d0 2 ) 6=3L « m 2 — 
K/i 2 , and M| c = {d 2 V/d6 2 ) e=Q « -m 2 + k/x 2 - Ca* 4 - Notice the sum rule: M 2 F + M 2 C = 
~( v 2qI j ' 2 ) 2 — 0, which guarantees that either or both of these curvatures must be negative 
(for all (j,). A third, mixed, phase (MX) — Zhang's 'spin-bag' phase — exists when both 
M 2 F and M 2 C are negative. In this case the minima are cos 26 MX = (M 2 C — M 2 F ) / (v2 Q — 
v 2s ) 2 , and the curvature at this minimum is M^ x = -2 M 2 C M 2 F /(M 2 C + M 2 F ). 

The critical doping which defines the boundary between these phases is found by 
solving the conditions M AF (fi AF ) = and M 2 c (/i 2 c± ) = 0, giving n\ F = n 2 sc _ + 0(e 2 ) = 
(m 2 j 'k) + 0(c 2 ), and n 2 sc+ = («/£) + 0(e). We see that the AF phase is predicted to extend 
only for dopings which are 0(e), i.e. a few %. The mixed phase, whose existence requires 
y? sc - < 1*afi occurs over a range of dopings which extend to 0(e 2 ) past the AF phase 
boundary, if it exists at all. Then one enters the SC phase, which persists over a wide 
0(1) range of dopings. These results break down for \i ~ /U sc +, since there the expansion 
in powers of fi 2 and 1/E 5 fail. 

For temperatures satisfying To(/i) < T < E$ SO (5) remains broken, but thermal 
transitions are possible between the SC and AF phases. An estimate for T (n) is obtained 
by finding the height of the potential barrier in V as a function of \i. For /j, chosen so that 
it is the AF phase which minimizes V, we find: T (AF) oc Harrier = \(M 2 F — M 2 C ) 
\(vp? — K[i 2 ). On the other hand, choosing \i so that V favours the SC phase gives: 
T (SC) OC V barri er = J(M| C - M 2 AF ). 

As discussed by Zhang, in each phase the four modes, 6> Q , group themselves into bona 
fide gapless Goldstone modes (having dispersion relation u 2 (k) = c^k 2 at low energies), 
and pseudo- Goldstone states (for which u 2 (k) = c 2 k 2 + e 2 ). The lagrangian given by (1) 
and (2) predicts the observables cq, c p and e for each phase. For brevity we record here 
only the lowest-order predictions for the AF and SC phases. 

In the AF phase two modes are gapless (magnons), with the other two forming an 
electrically charged pseudo-Goldstone doublet with respect to the unbroken SO (2). We 
find (co)af = (c p ) AF = c s = f s / ft + 0(e, fx 2 ). The gap is given as a function of doping by 
e 2 AF (^ 2 ) = (m 2 -^ 2 )/f 2 . 
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In the SC phase there is only one gapless mode, which is 'eaten' by the photon via 
the Anderson-Higgs mechanism, with the remaining three states forming a spin-triplet 
pseudo-Goldstone state. It is this state which Zhang argues persuasively [6] has been seen 
in neutron scattering experiments [7] in the spin-flip channel. The predictions in this case 
are: (c ) sc = (c p ) sc = c s , and e 2 sc = (-m 2 + Kfi 2 - ^ 4 )// t 2 . 

The final observables for which we present predictions are the electric and magnetic 
stiffness in the SC phase. Coupling to an electromagnetic potential, and eliminating the 
SC phase's Goldstone mode by an appropriate gauge choice, we find the term: CJ = 
\cl 2 e Aq — ia 2 f A 2 , with a E = qeft and a M = qef s . a E and a M represent the medium's 
electric and magnetic stiffness, normalized so that the penetration depth of electric and 
magnetic fields are respectively given by a" 1 and a" 1 . 

By eliminating the parameters m 2 , k, £, ft and f s we obtain several parameter- 
independent relations among the observables described above. We summarize these here: 



4» = %^ 



t^AF 



2 2 



Jl ( , \ £ sc( maX ) ( ,2 ,2 \/r>. 2 , 2^ 



L(0) £ 2 sc (max) (5) 



^ AF A*max 

.2 _ ..2 , ^ 



Vaf = Vsc- )> 
(cq)af = (c p )af = (c p )sc = (co)i 



a 



E 



In these expressions e 2 AF (0) = m 2 / f 2 denotes the AF gap energy at zero doping, and 
£| c (max) = (k 2 /4£/ 2 ) denotes the maximum gap size in the SC phase, which occurs when 

These are robust predictions of the assumed symmetry-breaking pattern regardless 
of the nature of the microscopic physics which is believed to be responsible, making their 
experimental verification particularly interesting. They receive calculable corrections in 
powers of e, which are most easily computed using the effective lagrangian presented here 
[4]. More predictions are possible should the MX phase arise, since the properties of the 



three Goldstone and one pseudo-Goldstone states in this phase are predictable in terms of 
the same parameters used here. 

We close with a speculation concerning how the 50(5) model might account for the 
anomalous electrical resistivity of the normal phase of the high-T c materials, which depend 
linearly on temperature (when doped with holes and near optimal doping). Within the 
SO (5) model this resistivity arises within a 'normal' phase below the >SO(5)-breaking 
scale, but at temperatures higher than the barrier between the AF and SC phases. The 
characteristic new feature of this 'normal' phase is the existence there of the light pseudo- 
Goldstone states, some of which carry electric charge. No such states are present in the 
normal phase of a BCS superconductor. In the rest of this section we wish to give a (naive) 
argument that if these light states play a significant role as carriers of electric current, then 
their scattering might give rise to a resistivity which is linear in T. 

To see why this is so, we estimate the temperature dependence of the resistivity by 
computing the temperature-dependence of the two-body rate for scattering of the appro- 
priate charge carriers from various quasiparticles within the material: p(T) ~ (ncrv) ~ 
nII|A| 2 . Here n describes the density of 'target' quasiparticles from which the charge 
carriers scatter, while A describes the amplitude for this scattering and IT represents the 
phase space available for the final states. The average is over the thermal distribution of 
scatterers, and the initial distribution of charge carriers. This kind of estimate reproduces 
the low-temperature T 5 dependence when applied to electron-phonon scattering, and the 
T 2 dependence when applied to electron-electron scattering. 

The basic observation now comes in two parts. First, notice that the electrically- 
charged, pseudo-Goldstone states of Zhang's SO (5) model obey a 'relativistic' kinematics 
since their gap energy is smaller than the temperatures of interest. Second, unlike low- 
energy phonons, the pseudo-Goldstone states need not interact only through derivative 
couplings, and so their scattering amplitude, A, from other particles need not be sup- 
pressed by powers of T at low temperatures. Combining these observations gives the 
standard estimate for the scattering rate of a 'relativistic' particle, which follows purely on 
dimensional grounds: p oc nII|A| 2 oc T. In this, admittedly simplistic, picture a resistivity 
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which is linear in T is seen as the signature that an appreciable part of the electromagnetic 
current is being carried by the light pseudo-Goldstone charge carriers. This motivates bet- 
ter understanding of the transport properties of these systems in this temperature range. 

The results of our analysis are encouraging, with the effective theory naturally ac- 
counting for the features of the high-T c phase diagram. The predictions of eqs. (5) make 
an even more quantitative test of the S'0(5)-invariant picture. Clearly much more remains 
to be done with this low-energy lagrangian, most notably including a systematic discus- 
sion of their transport properties, and an analysis of the response of the pseudo-Goldstone 
modes to various probes. 
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